USTC-ICTS-12-17 
November 2012 



On the angular distributions of r — > Kgn v t decay 



Dao-Neng Gao^ and Xian-Fu Wang* 

Interdisciplinary Center for Theoretical Study, University of Science and Technology of China, 

Hefei, Anhui 230026 China 



Abstract 



Based on experimental data by Bell Collaboration, we present a phenomenological 
analysis of the angular distributions in t~ — > Ks^~v T decay. Our study shows that the 
angular analysis could lead to some interesting observables, and the future experimental 
investigation of these observables might be very helpful to reveal the nature of the scalar 
components of the decay. New physics contributions from two-Higgs-doublet model to 
this decay have also been examined. 
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1 Introduction 



The t lepton is the only known lepton massive enough to decay into hadrons and its hadronic 
decays provide a very useful laboratory to test the low energy dynamics of the standard model 
(SM) [HI2]- With the increased experimental sensitivities achieved already or in the future, 
some interesting limits on possible new physics contributions to the r decay amplitudes 
could be expected; for instance, in the decays — > Kstt i/ T , upper limits on CP violation 
parameter from multi-Higgs-doublet models have been obtained by Belle Collaboration [3]. 

The decay r — > Kiri> T , which has the largest branching ratio of all Cabibbo-suppressed 
decays, could be a powerful probe of the strange sector of the weak charged current [U [5j [6]. 
Early investigations have established that the dominant contribution to the decay spectrum 
is from the if* (892) meson [TJ, while the scalar or tensor contributions, which are expected 
theoretically [EJIE], are not excluded experimentally JUfTU]. Recently, a precise measurement 
of t~ Ks^Vt decay has been done by Belle Collaboration [11] with high-statistics data, 
which showed that the K*(892) alone is not sufficient to describe the Kstc invariant mass 
spectrum and other states such as scalar (^(800), -K"q(1430), ...) or more vector resonances 
(K*(U10), K*(1680), ...) have to be included. Thus three models K*(800) + K*(892) + 
K*(U10), K*(800)+K*(892)+K*(U30), and K* (800)+ K* (892) +K* (1680) were introduced 
to fit the data, and the best description of the decay spectrum is achieved in the first 
two models One can expect that these two models will lead to the different angular 

distribution behavior since they have different scalar and vector components. This angular 
analysis has not been done yet experimentally. 

Motivated by the above new experimental data, in this paper we present a phenomeno- 
logical analysis of the angular distributions in the r~ — > Ks^~v T decay. We hope that this 
angular analysis would reveal the difference from the above models, which could further 
increase our understanding of the scalar contributions in this decay. 

The paper is organized as follows. In section 2, we will discuss the decay distributions 
of r~ — > Ks^Vti and some interesting observables will be introduced. In Section 3, a 
phenomenological analysis of the angular distribution is carried out, and new physics contri- 
butions from two-Higgs-doublet model to this decay will be examined. Finally, we summarize 
our results in Section 4. 



2 Decay distributions 

In the SM, the general invariant amplitude of r _ — > Kgit~v T can be decomposed as a product 
of a leptonic current and a hadronic current [B] 

M = ^£M„7». (1) 

Here Gf is the Fermi coupling constant and 6c is the Cabibbo angle. The leptonic current 
is given by 

M^ = u(p UT )j„(l-K)u(p T ), (2) 
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and the hadronic current can be parameterized by two form factors as 



J» = (Ksir-\vfu\0) = F v (s)[g 
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where Fy is the vector form factor, corresponding to the J p = 1~ component of the strange 
weak charged currents, and Fs is the scalar one, corresponding to the J p = + component. 
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In the limit of SU(3) symmetry, m\ 
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the vector current is conserved and Fs is zero. 



Now the differentiate decay rate, in terms of s, the ifs7r invariant mass squared, and 8, the 
angle between the three-momentum of Kg and the three momentum of r~ in the ifs7r~ rest 
frame, can be written as 
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and the phase space is given by 



ttik) 2 < s < m 2 



-1 < cos# < 1. 



After integrating over cos 8 in eq. 
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It is generally believed that this decay spectrum is dominated by the vector contributions 
(if* (892) resonances) [7]. Recent experimental results from Bell Collaboration [TT] have 
shown that the scalar contribution is necessary to fit the data, although it provides only 
a small contribution to the decay rate. It is found in [11] that two different combinations 
if *(800) + if* (892) + if* (1410) and if *(800) + if* (892) + if* (1430) can both describe the 
spectrum very well. In order to further understand the scalar form factor, one may carry out 
the angular analysis. From eq. (jl]) and only focusing on the angular part, the distribution 
can be written as 

io + ii cos 8 + i 2 cos 2 8, 
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Here (s m ; n , s max ) denotes the range of the integration over s, which could be the full phase 
space shown in eq. (jSJ) or some kinematical cuts on s. Using above equations, one can further 

get 

g + Acos0 + -{i-2g)cos 2 e, (ii) 



r d cos 9 

with 



r = 2(/ + \h) (12) 
is the decay rate. The constant term in (fTTl) can be expressed as 

S = ^ t (13) 

and the linear term in cos 9 

A=Y (14) 

will give a vanishing contribution to the decay rate after integrating cos# in the full phase 
space. However, this term can induce an interesting observable, called as the forward- 
backward asymmetry. Eq. fTH|) gives the integrated and normalized asymmetry. One can 
also define the differential asymmetry as [121 03] 

d 2 T \ r° ( d 2 T \ 

d cos 9—1 — — d cos 9 



.... . ds d cos 91 J-i\dsdcos9. 

(s) = ~' \, „ N <w \, ~' (15) 

d cos 9 + - — d cos 9 
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Thus together with eq. (TjJ, we have 

Pis) 

AMs) = 2 7 2 r S • (16) 

Note that the forward-backward asymmetries are generated from the interference between the 
scalar part and vector part amplitudes of the decay. In the case that the scalar contribution is 
not large, study of these asymmetries may be very important for us to extract the information 
on the scalar form factor. 



3 Phenomenological analysis 

In order to evaluate the quantities Q, A, and A FB defined in the previous section, we need 
the information on the hadronic form factors Fy and F$- Theoretically, due to the non- 
perturbative Quantum Chromodynamics (QCD) effects at low energies, there is no model 
independent way at the present. Therefore, various methods have been employed to study 
these form factors of r — > K-kv t decays, such as meson dominance models [6j E] and chiral 
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lagrangian including the resonances [151 HSJ El 03] • Phenomenologically, as shown in [TTJ [6] , 
these form factors could be parameterized by resonances contributions, which read 
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IS the s-dependent total width of the resonance. I = 1(0) for Kn in the P(S')-wave state 
and T 0R is the width of the resonance at its peak. (3, x, K , an d 7 in eqs.(TT7l) and (TTHT) are 
parameters for the fractions of the resonances, which would be determined by fitting the 
data of the hadronic invariant mass distribution. 

Actually these parameters have been determined in Ref. [TT] for two models with different 
combinations of the resonances, both of which can provide a good description of the decay 
spectrum. For the model I: K* (800) + if* (892) + if* (1410), 



0.075 ± 0.006, arg(/3) = 1.44 ± 0.15, 



1.57 ±0.23. 



(19) 



For the model II: (800) + if* (892) + #£(1430), the parameters could have two solutions 
from the data fit, which read (hereafter we call them as Model II-l and Model II-2, respec- 
tively) 



Model II-l : k 
Model II-2 : k 



1.27 ±0.22, | 7 | = 0.954 ±0.081, arg( 7 ) = 0.62 ± 0.34, (20) 

2.28 ±0.47, |7| = 1.92 ±0.20, arg( 7 ) = 4.03 ± 0.09. (21) 



It is obvious that there is no contribution from if* (1680) in these two models, \ in eq. ffTTl) 
is set to zero. 

Now we can calculate the quantity Q and the forward-backward asymmetry A. The 
results have been shown in Tables 1 and 2, respectively. It is found that, when these two 
quantities are evaluated for the full phase space, all the models give the consistent results. 
Model II-l and II-2 are consistent for Q even if we calculate it for different cuts on s. However, 
in the large s region, s >(1.2 GeV) 2 , it is expected that one might distinguish Model I and 
Model II from the values of Q. The asymmetry A shown in Table 2 is more interesting since 
its values will show the difference from these three cases for s above 1 GeV 2 . The similar 
conclusion can also be obtained from the differential asymmetry A-ps, which has been plotted 
in Figure 1. Afb's are almost consistent for three cases when s is below 1 GeV 2 , they however 
behave differently for s above 1 GeV 2 , in particular, Ap^ from Model II-2 could change sign 
for large s. 
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[m K + m^f ~(0.8 GeV) 2 
(0.8 GeV) 2 ~ (1.0 GeV) 2 
(1.0 GeV) 2 ~(1.2 GeV) 2 
(1.2 GeV) 2 ~(1.4 GeV) 2 



[\A GeV)< 



(m K + m n ) 2 ~ m 2 



Model I 



Model II- 1 



0.326±0.021 
0.257±0.001 
0.333±0.005 
0.398±0.003 
0.440±0.002 
0.269±0.002 



0.320±0.021 
0.257±0.001 
0.340±0.007 
0.436±0.006 
0.479±0.002 
0.270±0.003 



Model II-2 



0.343±0.039 
0.257±0.002 
0.339±0.008 
0.437±0.009 
0.479±0.005 
0.271±0.004 



Table 1: The values of Q defined in eq. ( JIB"]) are evaluated for different cuts on s in the above 
models [ ffl9l . (120]) . and fl21j) ]. The last line is for the full phase space. 
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(m K + m^f ~(0.8 GeV) 2 
(0.8 GeV) 2 ~ (1.0 GeV) 2 
(1.0 GeV) 2 ~(1.2 GeV) 2 
(1.2 GeV) 2 ~(1.4 GeV) 2 
(1.4 GeV) 2 ~m 2 
(rriK + m n ) 2 ~ m 2 



Model I 



Model II- 1 



-0.229±0.006 
-0.107±0.015 
-0.332±0.039 
-0.378±0.044 
-0.297±0.037 
-0.133±0.018 



-0.300±0.030 
-0.108±0.016 
-0.182±0.073 
-0.292±0.116 
-0.481±0.044 
-0.125±0.021 



Model II-2 



-0.217±0.022 
-0.105±0.032 
-0.413±0.075 
-0.042±0.110 
0.471±0.019 
-0.123±0.037 



Table 2: The integrated and normalized forward-backward asymmetry A are evaluated for 
different cuts on s in the above models [ (fl9~j) . ff20l . and ff2T]) ]. The last line is for the full 
phase space. 



Early studies tell us that the decay r — > Knv T is dominated by vector X*(892) con- 
tributions. The above different models [ ([19]) . f )20|) . and (|2T|) ] from Belle experiment have 
parameterized different scalar contents, therefore it is not surprising that analysis of the 
angular distribution of the decay will reveal some difference. This means that the future 
angular analysis from the high statistics experiment might be useful to distinguish the above 
models, and would help us to further understand the scalar components of the decay. 

In Ref. [13], the differential forward-backward asymmetry Afb has been studied in the 
two-Higgs-doublet model (2HDM) of type II with large tan/3 [18J. In this model, the charged 
Higgs exchange 

rH ± °f ■ Q m s m r tan 2 

L =—=sm9 c 2 M 1 + l5)rs(l - 7 5 )tt, (22) 

V 2 m H ± 

can lead to a tree level contribution to r — > K-kv T) which thus gives a new contribution to 
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Figure 1: The differential forward-backward asymmetry Afb is plotted as the function of 
\/s. The solid line is for Model I, the dashed-dotted line is for Model II- 1, and the dashed 
line is for Model 1 1- 2. 

the scalar form factor F$, 

KTow 1 m 2 K tan 2 B m s , , 

Fa =F S + -=— \^ — — . 23 

3 y/2 M 2 H± m s + m u K ' 

In general, this new contribution will be strongly suppressed by the large charged Higgs mass 
Mh±i which however may be substantially compensated by large tan/3. In terms of Fg ew , 
together with eq. ffT6l) . one can evaluate the new contribution to Afb, which has been plotted 
in Figure 2. We take F s in eq. ([23]) from Model I, and tan/3 /M H ± = 0.4 GeV- 1 ^ [20]- 
In Figure 2, we include the uncertainty from the present experimental fit for ^4fb ; which 
shows that the uncertainty could obscure the new physics signal, thus new physics searches 
through Ap-Q might be interesting in the future, however, might be not very significant at 
present. 



4 Summary 

Belle Collaboration reported a measurement of the decay spectrum of r~ — > i^57r~z/ T , and 
found two models, K o *(800) + if* (892) + if* (1410) and if *(800) + if* (892) + if* (1430), can 
both describe the spectrum very well. Based on these data, we carry out a phenomenological 
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Figure 2: The differential forward-backward asymmetry Ap-^ is plotted as the function of 
y/s. The red line is from Fg cw , the blue line is for Model I and grey region denotes the 
uncertainty from the parameters in Model I. 



analysis of the angular distributions in this decay. We find that the cos ^-dependence of 
the normalized spectrum 1 /TdT/d cos 9, see eq. tfTTj) . offers good opportunities to test the 
models. Therefore, the future experimental study of the angular distribution, in particular, 
the analysis of the forward-backward asymmetries, may be very useful to distinguish the 
above models, which would be helpful to reveal the nature of the scalar form factor of the 
decay. Possible new physics contributions from 2HDM to the Ap^ have also been analyzed, 
we found that the present experimental uncertainty may obscure the new physics signal. 
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